B NMR study of pure and lightly carbon doped MgB2 superconductors 
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We report a ^'^B NMR line shape and spin-lattice relaxation rate {l/{TiT)) study of pure and 
lightly carbon doped MgB2-a;Ci: for a; = 0, 0.02, and 0.04, in the vortex state and in magnetic field of 
23.5 kOe. We show that while pure MgB2 exhibits the magnetic field distribution from superposition 
of the normal and the Abrikosov state, slight replacement of boron with carbon unveils the magnetic 
field distribution of the pure Abrikosov state. This indicates a considerable increase of 1/^2 with 
carbon doping with respect to pure MgB2. The spin-lattice relaxation rate l/(TiT) demonstrates 
clearly the presence of a coherence peak right below Tc in pure MgB2, followed by a typical BCS 
decrease on cooling. However, at temperatures lower than ~ lOK strong deviation from the BCS 
behavior is observed, probably from residual contribution of the vortex dynamics. In the carbon 
doped systems both the coherence peak and the BCS temperature dependence of l/{TiT) weaken, 
an effect attributed to the gradual shrinking of the a hole cylinders of the Fermi surface with electron 
doping. 

PACS numbers: 74.25.-q., 74.72.-b, 76.60.-k, 76.60.Es 
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I. INTRODUCTION 



The recent discovery of superconductivity in MgB^i 
has motivated extensive experimental and theoretical 
work, as this material has the highest critical temper- 
ature among superconducting binary compounds (Tc ~ 
40K). MgB2 is isostructural and isoelectronic with in- 
tercalated graphite, and has a simple AlB2-type crys- 
tal structure (space group P6/mmm) consisting of al- 
ternating close-packed Mg layers and graphite-like hon- 
eycomb B layers. Most theoretical and experimental 
work on MgB2 has been devoted to the understand- 
ing of the mechanism responsible for its unusual super- 
conductivity. By now it is widely accepted that MgB2 
is a non-conventional BCS-type superconductor. First- 
principle calculationsSi^ show that superconductivity 
resides in two groups of bands: (i) the two strongly su- 
perconducting a bands and (ii) the two weakly super- 
conducting TT bands. A large number of experiments, in 
particular studies with Raman spectroscopji^, scanning 
tunneling microscopjiS*2iSiS, muon spin resonanceiSiii*i^, 
point-contact spectroscopjii^ti^ji^^i^ and specific-heal^ 
measurements support the notion that MgB2 is a dou- 
ble gap superconductor. On the basis of this evi- 
dence new physics has been unveiled. For example, 
theoreticaliSiiS and experimental^"-^-'^ works have shown 
that the anisotropy of the upper critical field 7^ 
H^2 / ~ 6 differs substantially from the anisotropy of 
the London penetration depth 7^ — \^ /X^^ w 1.2, at low 
temperatures. It has been also suggested that MgB2 pos- 
seses a peculiar vortex structure^, where vortices in the 



TT band are characterized by absence of localized states 
in the core, very large vortex core size, and strong core 
overlap. 

Many of the efforts have been focused in the study 
of doped MgB2 by partially substituting Mg or B sites 
with electron or hole dopants. Such substituents mod- 
ify (i) the crystal structure, and (ii) the electron den- 
sity of states at the Fermi level N{Ep), and therefore 
are extremely important in analyzing the band struc- 
ture and the microscopic nature of superconductivity 
in MgB2. A systematic study of Ala;Mgi_a;B2 sam- 
ples with Mg replaced by Al has already been under- 
taken by several groups iSSiSiSiS^ In particular, ^'^Al and 
"B NMR measurements in Al^Mgi^^B^^^iS^, pro- 
vided strong evidence that the Fermi surface is made 
of hole- type cr-bonding 2D cylindrical sheets, which fill 
up upon Al doping and collapse at a; « 0.55, together 
with hole-type and electron-type 3D 7r-bonding tubular 
networks, in accordance with theoretical predictionsiSi 
These results suggest anisotropic pairing and multi- 
gap superconductivity for MgB2. An alternative ap- 
proach is to introduce dopants into the B sublat- 
tice. Several groups have attempted to replace B with 
carboni".i2i28,29,30,3i,32^33,34,35,36,37,38,39 r^he results re- 
ported in the literature indicate that the degree of car- 
bon doping depends on the procedure and the starting 
material used in the synthesis of MgB2_a;C2:. Similarily 
with Al2,Mgi_a;B2, the replacement of boron with car- 
bon leads again to doping with one electron, and a sig- 
nificant depression of j'^jSSiSSi^iiSS, although the effect is 
more pronounced in the case of MgB2_xCa;. At the same 
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time, x-ray measurements have shown that MgB2_2;C2: 
exhibits a highly anisotropic lattice contraction with the 
lattice parameter a decreasing almost linearly with in- 
creasing doping level22i22iSi^2i22i, while the c parameter 
is virtually unaffected by the doping level, indicating that 
carbon is exclusively substituted in the boron honeycomb 
layer without affecting the interlayer interactions. Be- 
sides, data taken from tunneling'*'^, specific heat^*^, point- 
contact^, and muon SR measurementsii*i^ reveal that 
despite the homogeneous and random nature of C substi- 
tution the two-gap behavior is preserved in yigQ2-x^x- 
This is in contrast to theoretical predictions, accord- 
ing to which enhanced quasiparticle interband scattering 
should lead to an effective one-gap behavior in the doped 
samples^ and suggests that the interband scattering is 
selectively suppressed in y[^2-xGx^ 

In order to elucidate further the complicated electronic 
structure of MgB2, we have performed a comparative ^^B 
NMR study on pure MgB2 and carbon doped yi^2-xGx 
samples for x = 0.0, 0.02, and 0.04. Our results indicate 
a considerable increase of H^2 ^'^id strong supression of 
both the coherence peak and the BCS temperature de- 
pendence of l/(TiT) with carbon doping. These results 
are in agreement with recent ac-susceptibility and torque 
measurements^Sii^, which show strong enhancement of 
jy°2) passage from the clean to the dirty limit regime, 
and shrinking of the cr band hole-type parts of the Fermi 
surface with carbon doping. 



II. SAMPLE PREPARATION AND 
CHARACTERIZATION 

Polycrystalline pure MgB2 and carbon doped 
MgB2_a;C2; samples for x — 0.02 and 0.04 were prepared 
by reaction of Mg, amorphous B and carbon black 
at 900K for 2 h, as described in refi^ii^^ Synchrotron 
X-ray diffraction patterns, Raman and ac-susceptibility 
measurements of samples from the same batch are 
presented elsewerei^iiSSi^SiSi^ "B NMR line shape 
measurements of the central transition (—1/2 1/2) 
were performed on a home-built spectrometer operating 
in external magnetic field Hq = 23.5 kOe. The spectra 
were obtained from the Fourier transform of half of the 
echo, following a typical 7r/2 — r — tt spin-echo pulse 
sequence.^** The irradiation frequency was the same in 
all experiments. The spectral bandwidth around the 
irradiation frequency was « 150 kHz, which is enough to 
cover adequately the NMR signals at low temperatures. 
The ^^B Ti of the central line was determined by apply- 
ing a saturation recovery technique, and fitting with the 
two exponential relaxation function that is appropriate 
for / = 3/2 nuclei, m(t) = (M(oo) - M{t))/M{oo) = 
Q.lexp{-t/Ti) + 0.9ea;p(-6t/ri)i^ 



III. "B NMR MEASUREMENTS 
A. ^^B NMR line shapes and the mixed state 

When a type - II superconductor is placed into an ex- 
ternal magnetic field Hd < Ho < Hc2, at temperatures 
lower than the temperature of the second critical field 
Tc2, the magnetic field penetrates partially in the form 
of thin filaments of fiux, and generates a vortex lattice 
with a characteristic magnetic field distribution f{B). 
This can be excellently monitored by using magnetic res- 
onance techniques, like muon spectroscopy (for a review 
see ref4&) and NMRMiiLi^, as the resonance frequency 
is directly proportional to the local magnetic field. In 
case of an ideal hexagonal vortex lattice, a characteris- 
tic frequency distribution is formed with a peak corre- 
sponding to the saddle point of the magnetic field dis- 
tribution, which is located midway between two vortices, 
and two steps at the maximum (-ffmax) and minimum 
(-ffmin) fields^2iiS, corresponding to the vortex core and 
the center of a vortex triangle, respectively. This dis- 
tribution depends on the two fundamental length scales 
of superconductivity: the magnetic penetration length A 
(the variance of the frequency distribution is inversely 
proportional to A), and the coherence length ^, which 
defines the " extension" of a vortexi^Si^2iS It is thus pos- 
sible from the NMR line shapes to acquire information 
about the topology and dynamics of the vortex lattice, 
as well as parameters that reflect the symmetry of the 
superconducting pairing. 

The ^^B NMR spectra of the investigated pure and 
carbon doped MgB2 powder samples, in a magnetic field 
of 23.5 kOe, are shown as a function of temperature in 
Figures and|31 In the case of pure MgB2, the line 
shapes in the normal state remain unchanged down to 
Tc2, while for T < Tc2 a new peak develops at the low 
frequency side. In a recent work, the direct compari- 
son of the ^^B NMR line shapes with dc-magnetic sus- 
ceptibility measurements has shown that the frequency 
position and the shape of the low frequency peak by de- 
creasing temperature, follows the development of the vor- 
tex latticei^ In view of this fact, the persistence of the 
normal state signal down to the lowest measured tem- 
perature T = 5K, has been explained as manifestation of 
the high anisotropy in the upper critical field Hc2^^^ 
Considering that in pure MgB2 the anisotropy ratio is 
equal to 7 ?a 6 at low temperatures22*SiiS, it is expected 
that in external magnetic field of 23.5 kOe {H^2 < 2.35 
Tesla) a part of the randomly oriented grains remains in 
the normal state (unshifted peak) , while the rest is in the 
vortex state. This gives rise to superposition of NMR sig- 
nals from the vortex and the normal states. The above 
argument becomes obvious if we consider that the up- 
per critical field in each crystallite varies with the angle 
between the external magnetic field and the c— axis as, 
Hc2{0) = H^2{^ + (7| - 1) cos^ 9y^/^. According to this 
formula, only crystallites with Hq < Hc2{0) would give 
the characteristic signal of a type-II superconductor in 
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FIG. 1; ^'^B NMR line shapes as a function of temperature for 
pure MgB2 under a magnetic field H — 23.5 kOe. For clarity, 
each spectrum is normalized to its maximum intensity. 



the vortex state. If for a part of the grains, Ho > Hc2{d), 
then the NMR Une shape would be the sum of spec- 
tra coming from crystallites in the normal and the vor- 
tex state. Consequently, for Hq < H^2j the NMR line 
shapes will reflect solely the magnetic field distribution 
of the pure vortex state, whereas for H^2 < Hq only 
the line shape of the normal state will be observed. At 
this point we must note that the polycrystaline samples 
show a lower H^2j i^i comparison to single crystals, where 
H^2 ^ 30 kOe has been observed. This fact may be ex- 
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FIG. 2: "B NMR line shapes of MgBi.ggCo, ,02 as a function 
of temperature, under a magnetic field H — 23.5 kOe. For 
clarity, each spectrum is normalized to its maximum intensity. 

plained as showing that polycrystalline samples are more 
clean than single crystals and therefore exhibit a higher 

f 42.43 

It is thus quite remarkable that in the carbon-doped 
systems the normal state signal component disappears 
below Tc2, as clearly observed in Figures [3 and |3| In par- 
ticular, in the case of MgB1.98C0.02, the line shape starts 
to broaden below w 26 — 27K, while shifting to lower 
frequencies. The disappearance of the normal state sig- 
nal component indicates an enhancement (in comparison 
with pure MgB2) of H^2 above 23.5 kOe, thus giving rise 
to the pure vortex state at low temperatures. However, in 
the case of MgB1.96C0.04, hne shapes in the mixed state 
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FIG. 3: ^^B NMR line shapes of MgB1.96C0.04 as a function 
of temperature, under a magnetic field H = 23.5 kOe. For 
clarity, each spectrum is normalized to its maximum intensity. 



appear to be quite more symmetric, and less shifted than 
the line shapes of the x = 0.02 sample, which is a signa- 
ture of the gradual deformation of the vortex lattice by 
increasing carbon doping. Similar increase of H^2 ^-^d re- 
duction of the anisotropy have been previously observed 
in the mixed superconducting state of lightly Al-doped 
MgB2 samplesi^ 

The disappearance of the normal state signal in the 
carbon doped samples enables the exact calculation of the 
variance a (second moment) of the frequency distribution 
in the vortex state as a function of temperature. This is 
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FIG. 4: '^^B NMR line width as a function of temperature for 
pure MgB2 (filled circles are for the normal state and filled 
squares for the vortex state), MgB1.9gC0.02 (open stars), and 
MgBi.96Co.o4 (open circles). Lines are theoretical vs. 
T curves from refi^, after appropriate scaling for the follow- 
ing cases: (i) conventional BCS superconductivity, (ii) double 
gaped superconductivity with the magnetic field parallel to 
the c-axis (in the dirty limit), or to the ab-plane (in the clean 
limit), and (iii) double gaped superconductivity with the mag- 
netic field parallel to the c-axis (in the clean limit). 



important as the second moment of the NMR frequency 
distribution, is formally related to the magnetic penetra- 
tion depth A according to formula a w 7„0.0608((/)o/A^), 
where 7„ is the nuclear gyromagnetic ratio for boron, 
and 4>o is the flux quantum (= 2.07 x ■cm'^)^^. We 

note that in powder compounds with strong anisotropy 
one should take into consideration that a contains con- 
tributions from both in-plane Xab and perpendicular Ac, 
which may vary differently with temperature. According 
to the Ginzburg-Landau theory, clean superconductors 
with arbitrary gap anisotropy are characterized by a sin- 
gle anisotropy parameter, i.e. £,ab/£,c — ^ab/^c- Under 
this aspect one would expect that for MgB2, 7a = 75 ~ 6. 
However, recent theoretica l -'^^1^^ and experimentalSSiSi 
studies have shown that due to the two gap character of 
superconductivity in MgB2, 7a is nearly isotropic (« 1) 
at low temperatures, while it increases to ~ 2.6 at Tc- It 
is thus expected that at low temperatures a would reflect 
a mean 1/A^. Here, we would like to stress that there are 
certain factors which complicate the interpretation of cr in 
terms of the magnetic penetration depth: (i) as pointed 
out by Brandli^, various kinds of random distortions de- 
grade significantly the magnetic field distribution in the 
vortex state and (ii) in pure MgB2 the high frequency 
tail of the line shapes overlaps with the normal state sig- 
nal component, while in the carbon doped samples the 



5 



tail is strongly supressed. The latter effect might be pro- 
duced by the establishment of a disorder-driven Bragg 
glass state, right below T^2j which partially smooths out 
the details of the magnetic field distribution42i^ These 
problems allow only a rough estimation of the penetra- 
tion depth. For example, the slightly shorter second mo- 
ment of pure MgB2 at low temperatures in comparison 
with the carbon doped MgB1.98C0.02 sample, as shown in 
Figure 01 may reflect the stronger pinning effect by the 
carbon dopant Despite the above mentioned shortcom- 
ings, we would like to point out some features, which ap- 
pear to support the two-gap superconductivity scenario 
for MgB2. First of all, we notice that by assuming a mean 
A « I2OOA at low temperatures, we obtain a « 0.11 
MHz, which is not very far from the measured ct — ctjv 
values for the three studied samples (crjv is cr for the nor- 
mal state NMR signals), as shown in Figure 2| What is 
remarkable is that both carbon doped samples exhibit an 
almost linear temperature dependence of a with a small 
hump at ~ 17 K, which follows relatively well the pre- 
dicted temperature dependence for 1/A^(T), as described 
in a recently proposed two-gap model^. This may be 
considered as indication that two gap superconductivity 
is moderately affected at field 23.5 kOe, in agreement 
with recent point-contact spectroscopy data which show 
that at low temperatures the smaller tt gap survives up 
to fields « 60 kO o^^i^^ . As can be seen in Figure 01 there 
is relatively good coincidence between theory and experi- 
ment for two cases: (a) with the magnetic field along the 
c-axis in the dirty limit (predicted \ah ~ 1057A), and 
(b) with the magnetic field along the afo-axis in the clean 
limit (predicted Ac « 392 A). Nevertheless, the latter case 
can be excluded as it corresponds to order of magnitude 
broader vortex state NMR signals than those experimen- 
tally observed. On the basis of the above ^^B NMR line 
shape analysis, it may be concluded that carbon doping 
increases drastically the coherence length ^af,, while driv- 
ing the systems from the clean to the dirty limit regime, 
in agreement with recent observations42ii^ 
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B. 



NMR spin-lattice relaxation 



The temperature dependence of ^^B l/Ti for the pure 
MgB2 and MgB1.96C0.04 samples is presented in FigureEl 
In the case of pure MgB2 for T > 30K, Ti was measured 
on the peak of the normal state signal, and for T < 30K 
on two different spectral positions: (i) on the peak (sad- 
dle point) of the signal component from the vortex state, 
and (ii) on the peak of the signal component from the 
normal state, by applying the two exponential analysis 
described by Kotegawa et aliSSiS 

The biexponential nature of the nuclear magnetization 
recovery curves in this latter case is clearly demonstrated 
in Figure El which exhibits ^^B NMR saturation recov- 
ery curves at T = 5K for the normal state signal compo- 
nent (upper panel) and the vortex state signal component 
(lower panel). While the saturation recovery curve from 
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FIG. 5: "B NMR (l/Ti) relaxation rates as a function of 
temperature for pure MgB2 and carbon doped MgBi.96Co.o4- 



the vortex state signal can be fitted by considering a sin- 
gle Ti y = 3.5 s, in the case of the normal state signal 
component, a nice fit is only obtained by considering two 
relaxation times: a long one Ti y = 3.5 s, which follows 
the temperature dependence of Ti in the vortex state, 
and a short one Ti_„ = 0.2 s. Similar results have been 
reported in refsi ^^i^^ Apparently, the short Ti^n compo- 
nent belongs to the normal state NMR signal component 
from grains in the normal state, and will be neglected 
in the following. However, the corresponding l/(Ti^„T) 
values below Tc are enhanced over the Korringa value of 
the normal phase, as also pointed out in refi^, a phe- 
nomenon which is not yet understood. 

The evolution of the relaxation rate 1/Ti as a func- 
tion of temperature in the mixed superconducting state 
is shown in the upper panel of Figure [51 It is observed 
that above Tc « 30K, (l/(TiT)) is constant, with a value 
w 0.0056 K^^s^^, following a Korringa temperature de- 
pendence. Right below Tc, (l/(TiT)) forms a peak cen- 
tered at T w 28K. This peak is obviously attributed 
to the coherence (Hebel-Slichter) peakS and indicates 
a fully gaped s-wave pairing state. This effect has been 
recently demonstrated by Kotegawa et al. in a^^B NMR 
study of lightly doped Mgi_2;Alj;B2 samples;-- The solid 
curve in the inset of Figure [31 is the T dependence of 
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FIG. 6; ^^B NMR saturation recovery curves of pure MgB2 
at T = 5K, for measuring Ti of the normal signal compo- 
nent (uper panel) and of the vortex signal component (lower 
panel) . The vortex state signal was nicely fitted with a single 
exponential analysis. Fit of the normal signal component was 
only possible by considering a double Ti,„ and Ti^v analysis, 
where Ti.^ was set equal to the Ti value of the vortex state 
signal component (the dot-line is the best single exponential 
fit). 

l/{TiT), by using the typical s-wave model, where a phe- 
nomenological energy broadening function is assumed to 
be of a rectangle type with a width 26 and a height 1/26 
as presented by Hebeli^ By considering a superconduct- 
ing gap size of 2A/fcBTc « 4.5 and 6/A{0) w 1/3 in 
agreement with Kotegawa at alj^, a nice fit is obtained 
down to 15 K; however, (l/(TiT)) appears to saturate at 
lower temperatures. The origin of this effect is not quite 
clear. However, the opening of the second gap as possible 
origin of the deviation from the typical BCS behaviour 
is excluded as it does not gives the almost horizontal low 
temperature slope of the (l/(TiT)) vs T curve. Devi- 
ation from the BCS behaviour at low temperatures has 
been already reported in both conventional isotropi(42iS 
and anisotropic type II superconductors. 52i£i This phe- 



nomenon has been attributed to contribution in the relax- 
ation mechanism from the vortex dynamic^S*^, which 
becomes predominant at low temperatures, where the 
BCS contribution to the relaxation rates is very small. 
On the basis of this argumentation (l/(Ti)) in the mixed 
superconducting state of MgB2 is expected to be given by 
l/Ti = l/Ti,„ + l/ri,,™, where l/Ti,, and 1/Ti 

,vm, is the 

contribution from quasiparticle scattering and the vortex 
motion respectively. Considering the vortex fluctuations, 
we have l/Ti^^™ = •Tc/(1-|-w|t^), where ujl is the 
Larmor frequency and hj_ is the fluctuating field in the 
perpendicular direction. In the case of (uj^Tc) 3> 1 and 
by considering that h±^ ^ T, the above formula becomes 
1/^1, I'm ~ I'^h^^l uj\tc ^ T/tc, and l/Ti shows at low 
temperatures the linear T-dependence as in Figure [3 

At the same time, even a small ammount of carbon 
doping like x — 0.04 leads to significant suppression of 
both the coherence peak and of the rapid fall of l/(riT) 
below Tc, as observed in the lower pannel of Figure |S1 
First of all we would like to notice that the opening 
of the superconducting gap suppresses 1/Ti, while the 
sudden increase ol N{Ep) at Tc acts in the reverse way 
(l/(Tir) oc N^{Ef)), so the coherence peak in l/(riT) 
may be considered as the result of the competition be- 
tween these two factors. Hence, any change in the gap 
strength or N{Ep), would affect directly the coherence 
peak and the slope of the l/(TiT) vs T curve below Tc- 
It is also notable that (i) according to the literature non- 
magnetic inpurity- like carbon - scattering smears out the 
gap anisotropy and enhances the coherence peak£^®'^'§^, 
and (ii) the two-gap behaviour is retained upon light car- 
bon doping liii^Si^i*^ It is thus straightforward to con- 
clude that the supression of the coherence peak is not 
caused by impurity scattering, but is rather consequence 
of the shrinking of the a hole bands with carbon, i.e. 
electron doping. 



IV. CONCLUSION 

In conclusion, we have shown that even a small substi- 
tution of boron with carbon in MgB2 leads to substantial 
increase of ^^'^ unveils the magnetic field distribu- 
tion in the vortex state, as reflected in the ^^B NMR fre- 
quency distribution. The temperature dependence of the 
second moment of the NMR signal a vs. T shows a clear 
deviation from the BCS theory of the penetration depth 
A and correlates nicely with a two gap model for A. The 
H^2 enhancement indicates drastic reduction of S^ab and 
passage from the clean to the dirty limit with light carbon 
doping. According to the spin- lattice relaxation measure- 
ments, in pure MgB2 a strong coherence peak is observed 
close to Tc, followed by rapid decrease in the l/(TiT) vs. 
T plot, which is indicative of the s-type nature of super- 
conductivity. In carbon-doped systems, both effects faint 
out, presumably because of gradual electron filling and 
shrinking of the a hole-type cylindrical Fermi surfaces. 
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